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SUMMARY 

This document assesses the magnitude of velocity mea- 
surement errors that may arise when performing stereo 
PIV with the cameras viewing through a thick, refractive 
window and where the calibration is performed in one 
plane only. The effect of the window is to introduce a 
refractive error that increases with window thickness and 
the camera angle of incidence. The calibration should be 
performed while viewing through the test section win- 
dow; otherwise, a potentially significant error may be 
introduced that affects each velocity component differ- 
ently. However, even when the calibration is performed 
correctly, another error may arise during the stereo re- 
construction if the perspective angle determined for each 
camera does not account for the displacement of the light 
rays as they refract through the thick window. Care 
should be exercised when applying a single-plane cali- 
bration since certain implicit assumptions may, in fact, 
require conditions that are extremely difficult to meet in 
a practical laboratory environment. It is suggested that 
the effort expended to ensure this accuracy may be better 
expended in performing a more lengthy volumetric cali- 
bration procedure, which does not rely upon the assump- 
tions implicit in the single-plane method and avoids the 
need for the perspective angle to be calculated. 

NOMENCLATURE 

a distance from object plane (along normal) 

to center of lens (mm) 
c length of object field along x (mm) 

di distance between lens and image plane, 
measured along lens axis (mm) 
d 0 distance between lens and object plane, 
measured along lens axis (mm) 

/ focal length of lens (mm) 


H object height (mm) 

h image height (mm) 

i angle of incident light ray (deg) 

/ distance between lens and object centers (mm) 

M magnification 

n index of refraction 

p arbitrary point along object plane (mm) 

q coordinate normal to lens plane (mm) 

R distance between center of lens and 

Scheimpflug coincident point (deg) 
r angle of refracted light ray (deg) 

t window thickness (mm) 

X , Y coordinates in image plane (pixels) 

x,y coordinates in object plane (mm) 

z coordinate normal to object plane (mm) 

a ? angle between image and lens planes (deg) 

8 object plane distance (mm) 

8 perspective error (mm) 

0 perspective angle in x-z plane (deg) 

0 0 angle between object and lens planes (deg) 

X coordinate in lens plane (mm) 

a locations near lens measured along x (mm) 

Superscripts 

' perspective measurement 

Subscripts 

c geometric center of object 

1 position on image plane and coincident 

with the lens axis 

im image 

L left position, or reference to an arbitrary 

perspective view 
N normal to window 

o position on object plane and coincident 

with the lens axis 
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obj object 

R right position, or reference to an arbitrary 

perspective view 
r refracted 


Abbreviations 

PIV particle image velocimetry 


INTRODUCTION 

Errors in PIV are often categorized as either random or 
bias errors. The random errors are often consequences of 
aspects of the experiment that are difficult to control and 
contribute to a general uncertainty in the result, while bias 
errors are systematic and lead to a degree of over or un- 
der estimation of the velocity. Bias errors are the result 
of the characteristics of the illumination, particle size, im- 
age recording equipment, and image analysis procedures. 
For example, a typical bias error in PIV results from the 
tracer particle that is used to seed the flow, since a heavy, 
large particle will not follow the flow as faithfully as a 
smaller, more nearly neutrally buoyant particle. The bias 
errors may be of a predictable magnitude, but their reduc- 
tion requires a careful approach to each aspect of the PIV 
experiment and great care with the data analysis. Much 
effort has been expended in the analysis of these bias er- 
rors, and, for example, the error may be analyzed by the- 
oretical analysis, conducting experiments under carefully 
controlled conditions or by generating synthetic PIV im- 
ages of accurately known flow fields and subjecting these 
to a PIV analysis. The reader is referred to Raffel et al. 
(ref. 1), Lawson and Wu (ref. 2), and Soloff et al. (ref. 3) 
for examples of many such investigations. 

Control of the bias errors and improvement of the overall 
accuracy of PIV depends upon having an appropriately 
seeded flow, the correct choice of hardware for the ap- 
plication, and the careful application of the image anal- 
ysis algorithms. Another significant source of error lies 
with the calibration procedure, however, and the overall 
accuracy of the calibration depends on the quality of the 
calibration equipment and procedure as well as on the im- 
age analysis applied to the calibration images. A rigorous 
calibration for stereoscopic PIV is even more important 
than for two-component, single-camera PIV, as the dif- 
ference in the two measured velocity fields is required, 
which therefore amplifies any errors. 

The basic starting point for a calibration is a carefully 
manufactured calibration graticule that can be aligned 
precisely with and moved relative to the PIV laser light 
sheet to an accuracy of a few microns. Bjorkquist (ref. 
4) performed an analysis of the effect of calibration-plate 


offset and rotation with respect to the light sheet using 
synthetic PIV images, while Raffel et al. (ref. 1) assessed 
calibration-plate alignment errors using back projection 
of particle image recordings. Both studies showed that 
significant calibration errors could easily be the result of 
alignment errors equivalent to the laser sheet thickness. 

Besides the physical alignment of the calibration plate 
with the laser light sheet, it is essential that the optical 
conditions for the calibration are identical to those of the 
actual PIV test. This poses significant challenges for PIV 
in flows where there are significant changes in the refrac- 
tive index of the medium within the flow field, but also 
underlines how important it is that even seemingly triv- 
ial tasks such as imaging through windows in the test- 
ing apparatus should be done correctly. Indeed, it has 
been shown that a reduction of the image distortions of a 
stereo PIV test in a thick liquid medium can be achieved 
by placing a liquid-filled prism at the liquid/air interface 
(Prasad and Jensen (ref. 5)). 

The current analysis considers wind tunnel PIV tests, 
where the PIV imaging is performed through wind tunnel 
windows. It is noted that calibration errors for this prob- 
lem have been addressed by Soloff et al. (ref. 3), who 
showed that a volumetric calibration could adequately ac- 
count for the distortion introduced by refraction through 
the window. Of concern here, however, is what errors 
may result from the use of a standard calibration pro- 
cedure, where the calibration graticule is imaged in the 
plane of the laser sheet only. This is a widely used cali- 
bration procedure, and therefore an assessment of the er- 
ror is important. It is shown that a non-uniform image 
distortion occurs across the field of view that is propor- 
tional to the window thickness, and it is then suggested 
that this error is manifested in the stereo reconstruction 
of the velocity field. For test facilities in which thick 
windows are present, the calibration procedure should be 
modified to reduce this effect by performing a volumetric 
calibration or by compensating for the refractive effect on 
the perspective angle. 

REFRACTION EFFECT ON PERCEIVED 
PARTICLE DISPLACEMENT 

Principle of Stereo PIV 

It is firstly helpful to cover the basic theory for stereo- 
scopic PIV, so that the importance of the camera view, or 
perspective angle, may be understood. Full details of the 
theory may be found in Raffel et al. (ref. 1), for exam- 
ple. Figure 1 shows the relative positions of the object 
plane and an off-axis camera, along with a definition of 
the coordinate system. Figure 2 shows a schematic dia- 
gram of light rays for an off-axis PIV test (e.g., for stereo 
PIV), in which only perspective effects are considered. A 


2 




Figure 1 : Coordinate system and camera position. 


particle is recorded at two positions, 1 and 2. For con- 
venience, location 1 is coincident with the center of the 
object plane. The imaging system records the projection 
Ax' of the physical displacement A s_ onto the object plane. 
The perspective error £ x is the difference between the ac- 
tual displacement in the object plane Ax and the apparent, 
projected displacement Ax' such that 


Ax' = Ax + z x 

(i) 

e x = Aztan0 

(2) 


where 0 is the perspective angle in the x, z plane, which 
is the angle between the normal to the object plane and 
the line of the chief ray to the lens. The perspective angle 
varies along the light sheet, but it is usually close to the 
angle between the lens axis and the normal to the light 
sheet. Therefore, if there is a sizeable perspective angle, 
then the perspective error will be significant and the dis- 
placements Ax and A z cannot be determined unless the 
flow is imaged from two separate views of known per- 
spective angle. This leads to the so-called stereoscopic 
PIV, where the actual displacements can be determined 
by a stereo reconstruction of the two velocity fields by 
measuring the projected particle displacements from two 
different camera views. Regarding the two different cam- 
era views, L and R , of the same particle displacement in 
the flow, then there are two sets of measurements as fol- 
lows: 


Ax' r = Ax + Aztan0# 

( 3 ) 

Ax' l = Ax + Aztan0£ 

( 4 ) 


Combining these two equations gives the actual particle 
displacements in the flow as 



Figure 2: Notation for analyzing perspective error. 


Strictly speaking, the perspective effect of an uninten- 
tional rotational displacement of the cameras in the y-z 
plane also has to be considered, but in most applications 
it is small. While the above two expressions provide the 
displacements in the x and z directions, the displacement 
in the y direction is simply 

(6) 

since there is no perspective effect in this direction. 

Analysis of the Refractive Effect 

Now the presence of a test section window will be taken 
into account. The refractive effect will be considered for 
two limiting cases: (a) the plane of the light sheet is paral- 
lel to the glass window; (b) the plane of the light sheet is 
orthogonal to the glass window. The latter configuration 
is commonly used in, for example, cross-flow measure- 
ments. 

Object plane parallel to window- Figure 3 shows how 
the presence of a glass window can affect the light rays 
traveling from the particles to the image plane (i.e., the 
camera sensor). These rays are refracted as they pass 
through the window, which is idealized as a parallel sided 
block, and so they emerge from the window in a direc- 
tion parallel to the incident rays. The back projections of 
the refracted rays onto the object plane show the physi- 
cal distance between the incident and refracted rays. This 
displacement is given by 

8 = t (tan i — tan r) (7) 


Ax 

Az 


Ax^ tan 0# — Ax^ tan 0^ 
tan0# — tan0£ 

Ax' r — Ax' l 


tan0# — tan0£ 


( 5 ) 
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where i and r are the incident and refraction angles of the 
ray and t is the thickness of the window. The refraction 
angle is related to the incident angle by Snell’s law 


Refractive Medium 



Figure 3: Effect of refraction on particle displacement. 


ft = sin// sin r (8) 

where n is the refractive index and typically n = 1.5 for 
glass. It is helpful to consider what angles of incidence 
might be expected in a typical PIV setup. A 20-cm-wide 
area viewed normally from a distance of 2 m will yield 
incident angles at the window of no more than about 3 ° 
at the extremities of the viewing area. Typical stereo PIV 
bisecting angles for the case when both cameras are view- 
ing the same side of the light sheet are in the region of 
35°, so the refractive displacements in the stereo case will 
be at least 16 times greater than for viewing with a nom- 
inally zero perspective. The refracted particle displace- 
ment in the object plane is then given by 

Ax' r = Ax' + 8 2 — Si (9) 


so that 


Ax' r = Ax' + t [(tan /2 — tan/i) — (tanr 2 — tanri)] (10) 

where the second term in the above expression represents 
the refractive contribution to the displacement. Referring 
to figure 4, the local incident angle at the window corre- 
sponding to a point at an ordinate p along the object is 
given by 


tan / = 


(/sin0 c — c) + p 
/cos0 c 


( 11 ) 


where / is the distance from the lens to the middle of the 
object and 2c is the length of the object (i.e., the field of 
view). The refraction angle is straightforward to calculate 
from Snell’s law. The incident and refracted angles for 
an in-plane particle displacement Ax (equivalent to A p) 
can also be calculated, and hence the difference between 
the refracted particle displacement and the actual particle 
displacement can be found. 



Figure 4: Variable definitions with the light sheet parallel 
to the window. 


Figure 5 shows the relative difference between the ac- 
tual particle displacement and the refracted particle dis- 
placement for an example stereo view across a range of 
perspective angles, assuming a field of view of 2c = 0.2 
m, a viewing distance of / = 2 m, an in-plane particle 
displacement of Ax = 2 mm, and a window thickness of 
/ = 38 mm. The difference (based on equations 8, 10, 
and 11) represents an error and is presented as the ratio 
of the refractive error (difference between refracted and 
actual displacements) to the actual displacement. The er- 
ror is the least at p/c = 0, where the incident angles to 
the window are smallest. 

It is easily shown that the error itself is not strongly de- 
pendent upon the particle displacement, and the analysis 
above indicates that the significant sources of the refrac- 
tive error are the perspective angle 0 and the thickness t 
of the window. Since the error is directly proportional to 
the window thickness t, its effect has not been shown on 
the plot. The error increases with distance along the ob- 
ject since the angle of incidence increases in the direction 
p , and the error is positive for the same reason. 

The effects of field of view and viewing distance are sim- 
ply to change the range of variation of the error across 
the field of view as the incidence angle at the window 
changes. In stereo PIV the light sheet is viewed from two 
directions, and the perspective angles for the two cam- 
eras are frequently set to be approximately equal in mag- 
nitude. In the case where the two cameras view the flow 
from the same side of the light sheet, the refractive errors 
at the extremities of the field of view will be different 
for the two cameras. Conversely, in the case where the 
cameras view the flow from opposite sides of the light 
sheet, the refractive errors for the two cameras will be 
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> o .5 Glass window n= 1.5, t=3 8 mm 
o Particle displacement = 2 mm 
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Figure 5 : Particle displacement error with the light sheet 
parallel to the window. 



Figure 6: Variable definitions with the light sheet orthog- 
onal to the window. 


approximately the same. The particular refraction error 
under discussion here will not, of course, appear if the 
calibration is performed through the window, and the re- 
sults shown on the figure show that if the calibration is 
not performed while imaging through the window a sig- 
nificant error may be the result. 

Object plane orthogonal to window- Figure 6 shows 
the plane of the light sheet as being orthogonal to the 
window, and how the refraction of the light ray affects 
the position of the particle in the plane of the light sheet. 
The coordinate system for calculating the angle of inci- 
dence at the glass window is also shown. The difference 
between this light- sheet orientation and the previous case 
is that the shift of the light ray in the direction of the nor- 
mal to the glass window is important (i.e., 8#, as indi- 
cated). This shift is related to the refractive shift 8 by the 
expression 


Syv = 8/ tan/ (12) 

and so 

= t( 1 — tan r/ tan/) (13) 

The angle of incidence of a light ray at the window in this 
case is given by 


tan / = 


/cos0 c 

(/sin0 c — c) + p 


(14) 


and the rest of the analysis is the same as the parallel- 
window case. Figure 7 shows the error associated with 
refraction for the case with the light sheet orthogonal to 
the glass window and for the same basic conditions as 


for the parallel- window case. The key distinction here 
is that the angle of incidence decreases as p increases, 
which means that the difference between the perceived 
and actual displacements is negative and its magnitude 
decreases with increasing p. The magnitude of the error 
is significantly less than for the parallel orientation, pri- 
marily because the refractive shift 8^ is almost uniform 
along the light sheet. 


REFRACTION EFFECT ON THE CALIBRATION 
PROCEDURE 

Recall from the analysis in the previous section that the 
actual particle displacements for stereo PIV are given by 


Ax 
Ay 
A z 


Ax' l tan 0# — Ax R tan 0^ 
tan0# — tan 0^ 

2 

Ax’ r + Ax’ l 
tan0# — tan0/. 


(15) 


The PIV analysis and calibration provide the projected 
displacements (Ax f L , A y' L ) and (Ax R , A y R ). The perspec- 
tive angles 0 l and 0/? are used for the stereo reconstruc- 
tion to find the three velocity components, and these must 
be known to high accuracy for each camera as otherwise 
they will introduce an error into the stereo reconstruction. 
They may be determined by either (a) inferring them from 
the optical magnification of the system and using the sys- 
tem geometry for the Scheimpflug condition, or (b) per- 
forming geometric measurement of the object, lens, and 
cameras positions in space. An analysis of the refractive 
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Figure 7: Particle displacement error with the light sheet Figure 8: Scheimpflug condition and coordinate system 
orthogonal to the window. relative to the lens. 


error on these two methods is provided in the following 
sections. 


Determining the Perspective Angle from Optical 
Magnification and the Scheimpflug Condition 

The Scheimpflug condition is met when the object, lens, 
and image planes are coincident at some point in space. 
This provides the basis for the geometry indicated in fig- 
ure 8. From the geometry, the perspective angle 0 at an 
arbitrary point p along the object plane can be determined 
from 


tan 0 = 


d 0 sinQ 0 -p 0 +p 
d 0 cos 0 O 


(16) 


It is the local perspective angle that should be used for 
the stereo reconstruction. The important parameters to 
be found are those relating to where the lens axis crosses 
the object plane (0 O , p Q , and d Q ). Following the analysis 
of Prasad and Jensen (ref. 5), it is convenient to place a 
coordinate system (k,q) at the center of the lens. The A 
coordinate defines the position along the lens plane, and 
the q coordinate defines the distance of a point away from 
the lens plane in the direction of the object. The param- 
eter d Q is the distance to the object plane from the lens 
plane measured along the lens axis. The coordinates of a 
point on the object plane are given by 


q 0 bj =d 0 - r k tan0 (17) 

where the obj suffix indicates a point on the object plane. 
With the system correctly arranged in the Scheimpflug 
condition, all points on the object plane are brought into 


focus on the image plane by the lens. If qi m is the q co- 
ordinate position of the image of the object on the image 
plane, then the local magnification is given by 


M=- 


dim 


tfobj 


(18) 


The negative sign in the above expression is necessary 
because qi m is negative. When the image is in focus, the 
q coordinates of the object and image are related to the 
focal length / of the lens according to 


11 1 

f q ob j qim 


(19) 


The above two equations can be combined to give 


1 -\-M 

q °bj- f M 


( 20 ) 


If the object with a known set of overall physical dimen- 
sions is imaged, and the dimensions of the image can be 
measured accurately, then the local magnification M for 
any position p along the object can be determined, and 
hence the object distance d Q bj can be calculated. The co- 
ordinate A can be determined after additional considera- 
tion of the Scheimpflug geometry, but if it is considered 
that in a PIV calibration the coordinate p along the object 
is known at a number of locations, then the perspective 
angle Q 0 may be determined from 


A q Q bj = A/?sin0 o (21) 

where A q 0 bj is the change in the object coordinate q Q bj 
for a known increase Ap in the object coordinate p along 
the object plane (recall figure 8). It is therefore possible 
to find Q 0 if the magnification can be determined, which is 
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Figure 9: Calibration grid notation with object coordinate 
system. 
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Figure 10: Image of calibration grid with image coordi- 
nate system. 


5. Calculate the object distances q Q bj for the left and 
right hand sides of the image (using equation 20), 
and hence determine the change in lens plane to ob- 
ject distance A q 0 bj between the two ends of the ob- 
ject, taking care to respect the sign of A p. 

6. Calculate 0 O (using equation 21). 


Figure 1 1 shows an image of a calibration grid. The di- 
mensions of the extreme left and right hand sides of the 
image of the calibration grid are indicated (the pixel size 
in this case is 6.45 pm ), along with the object distance 
between the left and right hand sides. The left hand mag- 
nification is Ml = 0.0964, and the right hand magnifica- 
tion is Mr = 0.0891, which for the lens focal length of 
/ = 135 mm gives (using equation 20) object distances 
of <lobj,L = 1535 mm and q a bj,R = 1650 mm respectively. 
According to equation 21, this gives the perspective angle 
as 0 O = 55.23°. The above procedure could of course be 
carried out for a larger number of points across the cali- 
bration grid, in which case a linear fit between A q 0 bj and 
A p for the set of calibration points would be obtained. To 
find the local perspective angle, the distance d 0 and the 
coordinate p 0 have to be found. To do this, the above 
analysis can be repeated to find the image plane coordi- 
nate qi m and the physical distance along the sensor, which 
provides the angle oq between the image plane and the 
lens plane from exactly the same theory as above but ap- 
plied to the image plane. The distance R between the lens 
center and the Scheimpflug coincident point can be deter- 
mined since 


possible in a digital system if the pixel size of the sensor is 
known, or in a film system by direct measurement of the 
film recording. With reference to figures 9 and 10, which 
show a simple calibration grid and its image respectively, 
this is done as follows: 

1 . Measure the object heights Hl and Hr at the left and 
right hand sides of the calibration grid. 


d 0 = 7? tan 0 O 

di = 7? tana*- (22) 


where di is the distance of the image plane from the lnes 
center, measured along the lens axis, and by definition 


1 _ 1 1 

l~Jo + di 


(23) 


2. Measure the horizontal distance A p between the left 
and right sides of the grid. 

3. Determine the pixel positions (X,7) of the images 
of the objects at the left and right hand sides of the 
camera sensor. The image heights A Yr and A Yr in 
pixels are required. Then convert these to the physi- 
cal image heights Hr and Hr by multiplying the pixel 
heights by the pixel size. 

4. The magnification is M = h/H, so the magnifica- 
tions at the left and right hand sides of the image 


This then permits d Q to be calculated and effectively lo- 
cates the origin of the (A, q) coordinate system relative to 
the object plane, thereby allowing p Q to be found. For 
the present example, and with reference to data shown 
in figure 11, the above theory gives q^L = -148.01 mm 
and qim.R = -147.03 mm (so A qi m = 0.987 mm), and for 
the calibration grid image width of 7.721 mm the remain- 
ing Scheimpflug angle is a, = 7.34°. Note that calculation 
precision is essential here, so if the reader wishes to check 
the calculations, the raw data on figure 1 1 should be used 
to calculate the object plane parameters again to preserve 
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Ax = 1197 px (equivalent to 7.721 mm) 


Aq 0 bj, r = Aprsind, 


(24) 


AY l = 897 px 
h L = 5.786 mm • 
H l = 60 mm , 


Ap = 140 mm 



Figure 11: Image of a typical calibration grid. 



Figure 12: Virtual coordinate system for lens and object. 


precision, rather than simply using the object plane val- 
ues given above. Thus R = 1141 mm giving d 0 = 1643 
mm and p Q - 112 mm. 

The above analysis procedure requires that the calibra- 
tion object size be known, the image dimensions be mea- 
sured accurately, and the Scheimpflug condition be sat- 
isfied, and that each of these are achieved with care. It 
now becomes clear how the refraction through the win- 
dow can affect the calculation of the perspective angle, 
as the presence of the refracting window has the effect of 
shifting the position of the object. Figure 12 illustrates 
this, where the light ray from point (k,q 0 bj) actually ar- 
rives at the lens as if it were from the shifted virtual lo- 
cation (k n q 0 bj,r) with the glass window removed. The 
virtual location is shifted relative to the actual location 
by the refraction. When finding the perspective angle the 
refracted object positions must be used, so strictly speak- 
ing, 


where A q Q bj,r is determined from the above described im- 
age analysis without the need for any refraction correc- 
tion. To find Ap r , the back projection of the actual object 
displacement A p has to be found, and this is achieved 
using the same analysis as for the particle displacement 
shift above. This could be done iteratively by initially 
setting A p r = A p, finding 0 O , applying a refractive cor- 
rection to find A p r , and repeating the calculation until a 
suitable convergence is achieved. By applying the back 
projection, the refractive error on the distance d 0 is au- 
tomatically accounted for. Local refractive corrections 
would have to be applied to find the local perspective an- 
gle, and the procedure for this would be similar to above. 


A typical calibration graticule consists of a grid of points 
(a dot pattern or a line pattern), and given that the angle 
0 O is to be determined from the distance between the cali- 
bration points, the effect of the refractive shift is the same 
as the effect of refraction on the particle displacement for 
PIV. A typical grid pattern pitch is 5 mm, and figure 13 
shows the refractive error in the quantity Ap = 5 mm for 
the window thickness of 38 mm for the parallel and or- 
thogonal orientations for perspective angles of 35 ° and 
55°. The choice of these two perspective angles is useful 
because in the parallel case the perspective angle of 35 ° 
results in a perspective angle of 55° for the orthogonal 
light sheet/window orientation if the window orientation 
only was changed to the orthogonal configuration, leav- 
ing unchanged the camera and laser sheet positions. As 
with the PIV displacement error, the calibration stretch is 
of significantly smaller magnitude in the orthogonal ori- 
entation case. The error of 1.2% corresponds to an an- 
gular error of about 0.7° from a basic perspective of 35°, 
and given that tan 0 is required in the stereo reconstruc- 
tion, the significance of this would be approximately 3% 
error in the out-of-plane velocity estimation. For the or- 
thogonal orientation the error is less than half that of the 
parallel case. 


Determining the Perspective Angle from the System 
Geometry 

Figure 14 shows the geometry of the one half of the stereo 
PIV system for finding the perspective angle in the par- 
allel configuration. The local perspective angle at an ar- 
bitrary point along the object is found by the same basic 
method as described in the previous section, except here 
00 and p 0 are to be measured directly. It is difficult to 
ascertain exactly where the lens axis lies, and it is eas- 
iest in practice to measure the dimensions to the center 
of the object. Therefore, it is assumed in the remainder 
of this analysis that the angle 0 C is to be determined by 
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Figure 13: Calibration pitch error due to refraction. 


direct measurement. The geometric center of the cali- 
bration object is straightforward to define, but the center 
of the lens is not, so the accuracy of the angle 0 C will 
improve for large imaging distances. However, with the 
window present the lens axis on the lens side of the win- 
dow is shifted relative to the object side, so to find the 
perspective angle the shift must be accounted for and so 

tan0 c = a r /<2 (25) 

For a window thickness of 38 mm and a central perspec- 
tive angle of 35°, the shift is 10 mm in the direction in- 
dicated. It is clearly advantageous in this case to have 
as large a viewing distance as possible. For the viewing 
distance of 2 m, the error caused by the refractive shift is 
0.2°. This is much smaller than for the previous case, but 
assumes that the geometry can be measured accurately in 
the first place. 


DISCUSSION OF THEORETICAL ANALYSIS 

The effect of the refracting window is to introduce a 
stretching across the field of view, and in the case of 
stereo PIV, where the perspective is nominally in one di- 
rection only, the stretching will affect only one compo- 
nent of the projected displacement. If the calibration is 
performed with the window in place, then this stretching 
is automatically accounted for in the PIV analysis and 
theoretically no additional error arises (see below for a 
discussion of the stereo reconstruction error). However, 
if the calibration is not done with the window in place, 
or if the window is somehow altered between calibra- 
tion and execution of the PIV experiment, then an error 


Figure 14: Projected and actual distances from lens to 
object. 


is introduced into the PIV analysis. This error is pro- 
portional to the window thickness and depends upon the 
nominal perspective angle, and in a severe enough case 
the magnitude of the error could exceed 2% of the actual 
projected displacement. For the case when the perspec- 
tive is in one plane only (single-axis Scheimpflug), only 
one component of the projected particle displacement is 
affected, and this means that derived quantities such as 
vorticity and divergence would be affected. The basic 
measurement error might seem to be small, but in PIV 
there are many sources of error that are compounded, and 
so the calibration should be performed in an optical en- 
vironment that is as close to the actual experimental as 
possible. 

It is clear from the analysis that if a “standard” calibra- 
tion using a calibration plate aligned with the PIV light 
sheet is performed with the cameras viewing through the 
window, then no systematic error in the projected particle 
displacements is expected to arise. However, the veloc- 
ity field has to be reconstructed from the set of projected 
particle displacements, and a source of error appears to 
arise from the manner in which the perspective angle is 
calculated. This error may be significant, and could be 
corrected for by accounting for the refraction effect when 
the perspective angle is calculated. However, given the 
prior experimental evidence of the performance of volu- 
metric calibration methods (Soloff et al. (ref. 3)), it is 
perhaps better to use such a calibration procedure to start 
with. 

The benefit of applying a refractive-error correction in 
the stereo reconstruction would arise in a retrospective 
analysis of an existing set of data, where the expense 
of repeating the experiment would not be worth the re- 
covery of a few percent of accuracy. In this respect 
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the evaluation of the perspective angle is a potentially 
significant source of error in the analysis of stereo PIV 
data anyway. The method described earlier for deter- 
mining the perspective angle from optical magnification 
and the Scheimpflug condition, where no measurement 
of the camera and calibration geometry is required, as- 
sumes that the Scheimpflug condition has been satisfied, 
and this requires that the image be judged to be correctly 
in focus across the field of view. 

In the analysis presented for determining the perspective 
angle from optical magnification and the Scheimpflug 
condition, the change in image distance Ag ?m was less 
than 1 mm across the image sensor distance of about 7 
mm, meaning that great accuracy is required in the calcu- 
lation of the angle oq. This implies that the Scheimpflug 
condition has to be satisfied with great precision. Fur- 
thermore, the locations of the calibration grid points in 
the image need to be measured with great accuracy in 
order that all the various geometric parameters in the cal- 
culation are computed accurately. For example, even as- 
suming that the Scheimpflug condition was precisely met 
for the case shown in figure 1 1 , if the vertical locations of 
the centers of the dots on the image are subject to an error 
of only one pixel for each side, the perspective angle may 
be in error by A0 O = 3°, equivalent to an error in tan0 o 
of 10% for 0 O = 55.2°. Equivalent pixel errors produce 
a smaller error in Q 0 as the perspective angle falls. The 
error in tan0 o is about 6%for00 = 35°, and the error in 
the w component of velocity is of this magnitude. 

In the context of the current study, the effect of the re- 
fractive error with the standard calibration technique is to 
compound an error that is potentially significant anyway, 
and it is proposed here on the basis of this analysis that 
all stereo PIV experiments be conducted using a volu- 
metric calibration, where the calibration is performed by 
imaging the calibration grid in a number of closely spaced 
(x,y) planes. Using this calibration method would act to 
emphasize the care required in calibration for stereo PIV. 
In this calibration method the stereo reconstruction does 
not require the use of the perspective angle (see Soloff 
et al. (ref. 3)), as the (u,v,w) flow field is built up from 
the stereo PIV data using the three-dimensional calibra- 
tion. This is intrinsically more robust than the standard 
calibration as it does not make any assumptions about the 
focussing arrangement and does not require measurement 
of the relative positions of the camera and light sheet. 

Indeed the studies of Lawson and Wu (ref. 2) and Soloff 
et al. (ref. 3) both employed the volumetric calibration 
technique, but used different camera focussing arrange- 
ments. Lawson and Wu (ref. 2) used a simple arrange- 
ment with no lens/camera offset and the camera tilted 
with respect to the calibration plane. Sharp focus was 
achieved using a high lens f- number. On the other hand, 
Soloff et al. (ref. 3) performed experiments with the 


camera sensor, lens, and object planes all parallel to each 
other, achieving field-of-view overlap by using a trans- 
lational lens shift. This ensures good focussing at low 
f-number, but the joint field of view is relatively small. 
Note that neither of these studies used the Scheimpflug 
focussing arrangement. The basic experimental approach 
applied in these two studies was otherwise similar, but 
for mathematical representation of the calibration geom- 
etry in the image coordinates, Lawson and Wu (ref. 2) 
used a bicubic spline interpolation while Soloff et al. (ref. 
3) used a least-squares polynomial, with cubic depen- 
dence on the in-plane displacements and a quadratic de- 
pendence on the out-of-plane displacement. In neither 
case was any assumption made about the camera, lens, 
and calibration grid setup. In addition, as the volumet- 
ric calibration method can account for any combination 
of (Ax, Ay, A z) particle displacement (within the limits of 
the calibration) then the technique ought to be less sen- 
sitive to a slight misalignment of the laser sheet and cal- 
ibration grid, which is known to be a significant source 
of error in stereo PIV (Bjorkquist (ref. 4), Raffel et al. 
(ref. 1)). Bjorkquist (ref. 4) deliberately applied rotation 
and alignment displacements in an experimental simula- 
tion of light sheet/calibration plane misalignment in PIV, 
but was able to correct for the target misalignment using 
data from the volumetric calibration, which would not be 
possible using the standard calibration. 

Given the significance of the perspective angle in the 
stereo reconstruction when the standard calibration is 
used, stereo PIV experiments in the case where the per- 
spective angle is large in both the (x, z) and (y, z) planes 
suffer an even larger source of potential error than in the 
case analyzed here, even without the presence of the re- 
fracting window, and this is because all three velocity 
components require stereo reconstruction. The volumet- 
ric calibration procedure used in this case would intro- 
duce no additional error when compared to the standard 
procedure. Although the analysis performed in support 
of the current study has indicated that the standard cali- 
bration procedure leads to errors, there is no suggestion 
that the technique should be abandoned altogether. It is 
merely reinforced that the various assumptions made in 
the analysis of the calibration images require that great 
care be taken in executing the calibration procedure, oth- 
erwise the stereo perspective angle may be in error. How- 
ever, perhaps this degree of care is best directed, instead, 
at performing the volumetric calibration. The analysis of 
the effect of the refracting window has therefore gained 
additional value in that it has highlighted how the as- 
sumptions used in the standard calibration can lead to a 
significant contribution to both the random and bias er- 
rors referred to earlier. 
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EXPERIMENT TO VERIFY THE REFRACTION 
EFFECT DUE TO A GLASS WINDOW 


The effect of the refracting window is to introduce a non- 
uniform distortion across the field of view, and the cali- 
bration will account for this if (and only if) sufficient data 
points are taken across the field of view to allow an ade- 
quate curve fit to the distortion. Provided the calibration 
was performed with the same optical path as the actual 
experiment, there is no reason to suppose that the pres- 
ence of a thick, refracting window will have an adverse 
effect on the PIV as far as determining the projected par- 
ticle displacements Ax' is concerned. The stereo recon- 
struction may well suffer, however, and therefore the ac- 
curacy of the three velocity components may be degraded 
by an incomplete treatment by a disregard for the pres- 
ence of the window. This is because the perspective angle 
0 plays a fundamental role in the stereo reconstruction as 
follows: 


Ax 
A z 


Ax' r tan — Ax' l tan 0# 
tan0£ — tan0# 

Ax'j — Ax' r 
tan0£ — tan0# 


(26) 


For the displacement Ay the perspective effect is often 
very small, so the following expression usually suffices: 


Ay = 


4vL + Ayg 
2 


(27) 


A set of experiments to investigate the importance of the 
refraction by the glass window is as follows. No laser 
is required, but some means of representing PIV parti- 
cles in the object plane must be used. This can be done 
using, for example, a photograph printed on a suitable 
medium of a pattern that will produce an image represen- 
tative of particles in an actual PIV experiment, and ef- 
fectively simulates an illuminated slice of the flow. This 
must be aligned precisely with the calibration target, as 
in an actual PIV experiment. 


1 . Set up the stereo PIV with both cameras in the same 
horizontal plane and with both cameras on the same 
side of the simulated sheet of particles to investigate 
the refractive effect with the light sheet parallel to 
the window. 

2. Without glass window: 

(a) Perform calibration. 

(b) Move synthetic particle pattern in precisely 
measured sets of (Ax, Ay, A z) steps and obtain 
images. It is essential that steps in the x, y, 
or z directions only are performed in addition 


to steps based on combinations of those direc- 
tions, as this will permit the effects of stereo re- 
construction to be isolated. Note that relatively 
few measurements of Ay displacement are re- 
quired, as in this arrangement the perspective 
effect is small. 

(c) Perform PIV analysis of data, and compare re- 
sults with actual displacements. The analysis 
of the PIV results should proceed carefully. 

3. Repeat procedure 2 with glass window. 

4. Repeat entire procedure from 2 with different stereo 
separation angles. 

5. Repeat entire procedure from 1 with a simulated par- 
ticle field and window in the orthogonal orientation. 

A data analysis to isolate calibration distortion from 
stereo reconstruction error would then need to consider 
at least the following points: 

1. For Ay displacements only (i.e., with other displace- 
ments zero), compare actual data with PIV data with 
and without the window. This will permit the theo- 
retically small refraction distortion in the y direction 
to be assessed. 

2. For A z displacements with a constant Ax displace- 
ment, compare A z displacements from the PIV anal- 
ysis with the actual values with and without the win- 
dow. The Ax' values (i.e., PIV analysis data be- 
fore the stereo reconstruction has been performed) 
should also be taken, so that the calibration distor- 
tion error can be isolated from the stereo reconstruc- 
tion error (i.e., the contribution due to error in per- 
spective angle 0). 

3. For Ax displacements only (i.e., other displacements 
zero) compare Ax' displacements from PIV with the 
actual Ax displacement with and without the win- 
dow. This will allow the effect of the perspective 
angle (and hence a larger refractive effect) to be seen 
in the refractive distortion on the x direction. 

An alternative to translating the simulated particle field 
by Ay would be to rotate the field about its z axis instead. 
This would simulate a rotational flow field with uniform 
vorticity (zero vorticity in this case) and would show how 
the process is affected when there are displacement gra- 
dients in the x and y directions. For the experiment just 
described, and given that only the Ax displacement would 
be significantly affected by the presence of the window, 
any error in the calculated vorticity could be attributed 
to a displacement gradient that is analogous to the du/dy 
contribution to vorticity. In other words, a truly rotational 
motion of the simulated particle field would result in an 
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erroneous vorticity with a magnitude that depends on the 
refractive error caused by the window. 


CONCLUSIONS 

An analysis of the error in stereo PIV data due to imag- 
ing through a glass window has been performed. When a 
standard calibration based on aligning a calibration plate 
with the PIV light sheet is performed, the error in the 
PIV analysis arises when the stereo reconstruction is per- 
formed as a result of an incorrect calculation or measure- 
ment of the perspective angle. It is indicated how the 
stereo reconstruction error may be isolated using an ideal- 
ized stereo PIV experiment. The perspective angle mea- 
surement may be significantly in error anyway, and it is 
suggested that the volumetric calibration method be used 
to avoid this source of systematic error. 
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